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Abstract 

A  new  equation  of  state  (EOS)  is  established  to  describe  water  activities  of  polymeric  membrane/water  systems  in  fuel  cell.  It  is  developed 
based  on  the  modified  perturbed  hard- sphere-chain  (PHSC)  EOS  by  introducing  a  new  perturbation  equation  that  is  obtained  from  the  generalized 
Lennard-Jones  (GLJ)  potential  function  based  on  a  statistical-mechanical  relationship.  Experimentally  observed  water  activities  of  polymeric 
electrolyte/water  systems  are  interpreted  by  the  proposed  model.  The  values  calculated  from  the  proposed  model  are  in  a  good  agreement  with  the 
experimental  data  for  given  systems. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  fuel  cell,  first  demonstrated  by  Sir  William  Grove  in 
1839,  has  taken  longer  than  most,  despite  the  promise  of  clean 
and  efficient  power  generation.  Since  DuPont  showed  that  per- 
fluorosulfonic  acid  ion-exchange  membranes  have  significant 
ionic  conductivities  and  can  be  use  in  the  fuel  cell,  considerable 
research  has  been  directed  toward  the  development  of  polymer 
electrolytes  having  high  ion  conductivity  as  well  as  mechanical, 
thermal,  and  chemical  stability.  One  important  factor  that  can 
affect  the  performance  of  fuel  cell  is  the  hydration  of  perfluoro- 
sulfonic  acid  membranes  [1-5].  Membranes  can  neither  be  too 
dry,  in  which  case  the  proton  conductivity  decreases  nor  too  wet, 
in  which  case  electrode  flooding  may  result.  To  solve  the  above- 
mentioned  problem  of  solid  polymer  electrolyte  (SPE)  for  fuel 
cell  application,  a  general  thermodynamic  model  to  be  able  to 
predict  phase  equilibria  for  these  systems  reliably  is  required. 
During  the  last  half-century,  there  are  various  theoretical  and 
experimental  approaches  to  investigate  the  structures  and  ther¬ 
modynamic  properties  of  polymers  and  their  solvents. 
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EOS  is  very  efficient  tool  to  understand  properties  of  polymer 
solutions.  Molecular-based  thermodynamic  models  for  describ¬ 
ing  vapor-liquid  equilibria  (VLE)  can  be  divided  into  four  cate¬ 
gories,  each  corresponding  to  a  particular  statistical  mechanical 
framework:  incompressible-lattice  models,  generalized  van  der 
Waals  partition-function  theories,  compressible-lattice  models, 
and  off-lattice  (continuous- space)  models  of  chain  fluids. 

Among  several  kinds  of  equations,  the  cubic  equation  of  state 
is  the  simplest  type,  which  is  capable  of  representing  both  liq¬ 
uid  and  vapor  phases.  Unfortunately,  simple  cubic  EOS  is  only 
applicable  to  the  calculation  of  fluid  phase  equilibria  in  sim¬ 
ple  systems.  Dickman  and  Hall  [6]  were  the  first  to  present 
an  EOS  for  hard-sphere  chains.  They  obtained  the  pressure 
through  the  formulation  of  chain-segment  insertion  probabili¬ 
ties  in  a  manner  analogous  to  that  used  to  derive  the  lattice- 
based  Flory-Huggins  theory.  Later,  Honnell  and  Hall  [7]  derived 
a  new  EOS  utilizing  the  Carnahan- Starling  equation  [8]  and 
the  Tildesley-Streett  [9]  dimer  equation.  These  equations  have 
been  tested  against  computer- simulation  results  for  hard- sphere 
chains  [10]  and  their  binary  mixtures  [11,12]  also  extended 
the  above  equation  to  include  square-well  chains  where  non- 
bonded  segments  interact  through  a  square- well  intermolecular 
potential. 

In  general,  hard-sphere  chains  can  provide  a  useful  reference 
system  for  representing  a  simple  model  of  polymers  [13].  Chiew 
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solved  the  Percus-Yevick  integral  EOS  coupled  with  chain  con¬ 
nectivity  for  mixtures  of  hard- sphere  chains  analytically.  How¬ 
ever,  as  Chiew  used  different  radial  distribution  functions  in 
bonding  term  and  non-bonding  term  [13],  Chiew’ s  equation  did 
not  show  the  relation  between  the  bonding  and  non-bonding 
terms. 

Song  et  al.  [13]  proposed  PHSC  EOS  for  the  purpose  of 
calculating  phase  equilibria  in  solutions  including  solvents  and 
polymers  by  generalizing  Chiew ’s  equation  through  reformu¬ 
lating  the  radial  distribution  functions  of  hard-sphere  mix¬ 
tures  at  contact  in  both  the  bonding  and  non-bonding  terms. 
It  gives  us  very  useful  information  in  investigating  properties 
of  fluids,  including  mixtures  of  chain-like  molecules.  It  con¬ 
tains  all  the  information  about  substances  and  any  equilibrium 
property  can  be  obtained  from  any  one  of  the  fundamental 
equations  of  state  by  applying  the  appropriate  thermodynamic 
relations. 

This  model  has  three  parameters  characterizing  normal  fluids 
including  polymers:  the  number  of  single  hard  spheres  per  chain 
molecule  (segment  number)  r,  the  separation  distance  between 
segment  centers  (segment  size)  a ,  and  the  segment-pair  inter¬ 
action  energy  elk. 

In  this  work,  we  employed  the  reference  equation  for  hard 
spheres  proposed  by  Hu  et  al.  [14].  To  describe  the  attrac¬ 
tive  forces  between  hard  spheres,  the  square-well  (SW)  fluid 
model  has  been  extensively  studied  until  recently  because  it 
is  useful  as  the  first  approximation  to  the  understanding  of 
properties  of  real  fluids  with  spherically  shaped  molecules.  In 
fact,  properties  of  real  fluids  can  be  approximated  by  those 
of  an  SW  fluid  (SWF)  with  suitably  chosen  potential  depth 
and  range.  Moreover,  thermodynamic  properties  of  an  SWF 
are  easier  to  determine  theoretically  than  those  of  other  flu¬ 
ids  having  more  complicated  potentials.  This,  in  particular, 
is  the  case  for  perturbation  theories,  from  which  thermody¬ 
namic  properties  of  SWFs  can  be  obtained  analytically  [15-17]. 
However,  the  SW  model  appears  to  be  a  rather  crude  approx¬ 
imation  of  a  real  fluid.  To  obtain  thermodynamic  properties 
of  a  real  fluid  from  those  of  a  fluid  with  a  model  poten¬ 
tial,  we  used  the  GLJ  potential  model  established  by  Kim  and 
Bae  [18]. 

Thermodynamic  properties  investigated  in  this  study  are 
water  activities  of  SPE/water  systems.  We  develop  a  new  ther¬ 
modynamic  framework  based  on  a  MPHSC  model  wherein 
proposed  by  Hu  et  al.  [14]  and  Kim  and  Bae  [18]  into  the  refer¬ 
ence  term  and  accounts  for  the  attractive  force  as  a  perturbation 
term.  We  also  calculate  the  water  activity  by  proposed  model 
and  GLJ  potential  in  PEM  fuel  cell. 

2.  Model  development 

2.7.  PHSC  EOS 

Derivation  of  EOS  for  mixtures  follows  a  rigorous  first  order 
statistical-mechanical  perturbation  theory  based  on  a  mixture 
of  hard-sphere  chains  as  the  reference  system,  and  closely  par¬ 
allels  that  for  pure  fluids.  Only  the  main  steps  need  to  be 
presented. 


The  van  der  Waals-type  equation  of  state  has  the  general  form 
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In  this  work,  we  employ  the  reference  equation  proposed  by 
Hu  et  al.  [14] 
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where  P  is  the  pressure,  p  =  N/V  the  number  density  (N  is  the 
number  of  molecules  and  V  is  the  volume),  and  k  is  the  Boltz¬ 
mann  constant  and  a ,  /3,  y  parameters  are  as  follows: 
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(72  =  0.45696,  Z?2  =  2.10386,  c2  =  1.75503,  a3  =  -0.74745, 

7>3  =  3.49695,  6*3  =  4.83207,  where  rj  is  defined  as  the  packing 
fraction  given  by 
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where  b  reflect  the  second  virial  coefficient  of  hard  spheres 
2j x  o 
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In  Eq.  (1),  the  perturbation  term,  which  accounts  for  the 
attractive  force,  can  be  generally  expressed  as  below. 
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where  U  is  the  perturbation  energy  that  can  be  obtained  from  the 
following  relationship: 
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where  W(r)  is  the  pair  potential  function. 

The  GLJ  potential  function  has  the  form  of  [18], 
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where  r  is  the  center-to-center  distance  between  two  adjacent 
molecules,  a  the  diameter,  e  the  minimum  potential  energy,  T  = 
kT/e  the  reduced  temperature  with  Boltzman  constant  k,  and 
A(T)  is  the  temperature-dependent  parameter  that  is  determined 
from  the  computer-simulation  data  for  the  compressibility  factor. 
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Combining  Eq.  (8)  with  Eq.  (9)  yields  the  perturbation  equa¬ 
tion  as  follows: 
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Then  the  final  form  of  EOS  can  be  written  as 
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The  EOS  for  the  mixture  is  obtained  by  extending  Eq.  (12) 
directly 
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where  x/  =  NJN  is  the  mole  fraction  of  molecule  i  and  r,  is  the 


number  of  hard  spheres  and  B ,  E ,  F  are  given  by, 


To  transform  the  number  basis  to  segment  basis  equations, 
more  useful  variables  should  be  introduced  for  polymer  mix¬ 
tures.  These  variables  are  segment  density,  pr  =  Nr/V ,  and  seg¬ 
ment  fraction  0;  =  Af;r0/Vr,  where  Nr  is  the  total  number  of 
segments  in  the  system  with  volume  V. 


m 


Nr  =  XV 


Ti 


(19) 


pkT 


In  these  terms,  the  final  form  of  the  new  equation  of  state  is 
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2.2.  Thermodynamic  functions 

The  Helmholtz  energy  can  be  calculated  from  the  pressure- 
explicit  EOS  [19] 
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The  additivity  of  hard- sphere  diameters  allows  the  expression 


for  bij  as  follows: 
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The  pair  potential  between  unlike  segments  Sy  and  size 
parameter  o y  can  be  determined  by 
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2.3.  Chemical  potentials 

The  chemical  potential  is  defined  as: 
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(a) 
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CF, 


(b) 


Fig.  1 .  Repeat  unit  structures  of  perfluorinated  proton  exchange  membrane  with 
sulfonic  acid  functional  groups,  (a)  Nation  (x  =  5-1 1)  and  (b)  Flemion  (m-  0  or 
1,  n-  2-5,  x=  1.5-14). 


and  activity  can  be  obtained  as 


Table  1 


Lists  of  characteristic  parameters  for  the  proposed  model 


Polymer 

Characteristic  parameters 

r 

a 

s/k 

Nation  117/water 

2735.2 

9.8 

1540 

Nation  125 /water 

2682 

10.82 

1585 

Flemion  4/water 

4326 

10.20 

1610 

Flemion  12/water 

2745 

8.84 

1442 

Aciflex  4/water 

2725 

10.40 

2400 

tively.  The  values  of  the  three  characteristic  parameters  are  listed 
in  Table  1 . 

All  results  show  that  the  water  activities  increase  as  the  equiv¬ 
alent  weight  of  the  polymer  decreases,  or  in  other  words,  as  the 
number  of  exchange  sites  increases.  The  calculated  curves  for 
the  different  polymer  electrolytes  show  the  same  trends  because 
all  the  commercial  SPEs  have  similar  chemical  structures  (as 
shown  in  Fig.  1). 
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3.  Results  and  discussion 

We  developed  a  new  thermodynamic  frame  work  based  on 
combining  the  reference  model  by  Hu  et  al.  [14]  and  the  pertur¬ 
bation  model  by  Kim  and  Bae  [18]  to  describe  water  activities 
of  perfluorosulfonic  acid  polymer  electrolyte  systems — Nation 
and  Flemion.  Their  structures  are  illustrated  in  Fig.  F 

The  model  has  three  characteristic  parameters:  the  num¬ 
ber  of  single  hard  spheres  per  chain  molecule  r,  the  segment 
size  a  and  the  non-bonded  segment-pair  interaction  energy 
s.  For  i  component  (water)  and  j  component  (polymer),  the 
pure  component  parameters  for  polymers  (rj,  oj ,  sj)  and  mix¬ 
ture  parameters  (r^,  (x,y,  Sij)  are  obtained  by  fitting  simulation 
results. 

In  Fig.  2,  the  solid  and  dotted  lines  are  calculated  from  the 
proposed  model  and  MPHSC  equation,  respectively  [20] .  In  the 
entire  density  region,  MPHSC  EOS  shows  slight  deviations  from 
the  experimental  data,  however  our  model  shows  good  agree¬ 
ment  with  the  data.  The  obtained  characteristic  parameters  for 
water  are  V water  =  F75,  O water  =  2.81,  £water/^=  680.2. 

Figs.  3a,  b,  4a,  b,  and  5a  represent  calculated  water  activities 
of  Nation  117,  125  (DuPont,  EW=  1100,  1200),  Flemion  4,  12 
(Asahai  Glass,  EW  =  890,  920),  and  Aciflex  4  (Asahai  Chemi¬ 
cal,  EW  =  980),  respectively.  Dark  squares  are  experimental  data 
reported  by  Hinatsu  et  al.  [21].  The  solid  and  dotted  lines  are 
calculated  by  the  proposed  model  and  MPHSC  equation,  respec- 


Fig.  2.  Vapor/liquid  equilibrium  experimental  data  for  water.  The  solid  line  and 
dotted  line  are  calculated  by  this  work  and  the  modified  PHSC.  Open  squares 
are  experimental  data  [20] . 
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Fig.  3.  Water  activities  for  the  Nation  (E.I.  DuPont  de  Nemours  and  Company )/water  system  at  353.15  K.  (a)  Nafin  1 17  (EW  =  1 100)  and  (b)  Nation  125  (EW  =  1200). 
The  solid  line  is  calculated  from  this  work.  The  dotted  line  is  calculated  from  modified  PHSC.  Dark  squares  are  experimental  data  [21]. 


0.45  0.50  0.55  0.60  0.65  0.70  0.75  0.80  0.85 


(b)  Weight  fraction  of  Water 


Fig.  4.  Water  activities  for  the  Flemion  (Asahai  Glass)/water  system  at  353.15  K:  (a)  Flemion  4  (EW  =  890)  and  (b)  Flemion  12  (EW  =  920).  The  solid  line  is 
calculated  from  this  work.  The  dotted  line  is  calculated  from  modified  PHSC.  Dark  squares  are  experimental  data  [21]. 


on  a  statistical-mechanical  relationship.  The  proposed  model 
includes  the  short-range  interaction  contribution  and  the  long- 
range  interaction  contribution  which  depends  on  temperature. 
The  characteristic  parameters  are  calculated  from  the  experi¬ 
mental  data  so  that  we  are  able  to  predict  the  water  activities  of 
fluorinated  solid  polymer  electrolyte/water  systems. 
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Fig.  5.  Water  activities  for  the  Aciplex-4  (EW  =  980)  (Asahi  Chemical  Industry 
Co.)/water  system  at  353.15  K.  The  solid  line  is  calculated  from  this  work.  The 
dotted  line  is  calculated  from  modified  PHSC.  Dark  squares  are  experimental 
data  [21]. 

4.  Conclusion 

In  this  study,  we  propose  a  new  thermodynamic  model 
based  on  the  MPHSC  for  SPE/water  mixtures  by  introduc¬ 
ing  a  new  perturbation  equation  that  is  obtained  from  the 
generalized  Lennard- Jones  (GLJ)  potential  function  based 
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